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Abstract 



An detailed study of the thermodynamical state of nuclear matter in 
^ ■ transport calculations of heavy-ion reactions is presented. In particu- 

'nJ" , lar we determine temperatures from an analysis of the local momentum 

space distribution on one hand, and from a fit to fragment energy spec- 
fSj ' tra in terms of a blast model with radial flow and temperature on the 

O . other. We apply this to spectator and participant matter. In spectator 

we find regions of spinodal instability with temperatures and densities 
which are consistent with experiments. In the participant we find dif- 
ferent temperatures for different fragment masses, indicating that the 
fragments are not emitted from a source in thermal and chemical equi- 



^ ' librium. 



1 Introduction 



H I One of the challenges in the study of heavy-ion collisions is the understand- 

ing of multifragmentation in relation to liquid-gas phase transitions. In spite 
of many experimental and theoretical efforts these processes have not been 
fully understood yet. This is largely due to the fact that a heavy ion colli- 
sion is a dynamical process where the state of nuclear matter varies strongly 
in space and time, and which during much of the reaction is not in global 
or even local equilibrium. This is, e.g. seen by looking at local momentum 
space distributions in transport calculations [|l], H which are found to be highly 
anisotropic even during the compression phase of the collision. Only at the 
later stages of the reaction the local momentum distributions become more and 
more thermalized without neccessarily leading to global thermal equilibrium. 
Therefore, non-equilibrium effects are important for a reliable description of 



heavy ion colhsions |1|, ^. The influence of these non-equihbrium effects on 
the determination of the equation of state of nuclear matter has been dis- 
cussed in Ref . [^ . In this contribution we concentrate on the question of the 
applicability of thermodynamical concepts in the non-equilibrium situation of 
heavy ion collisions [Q and, in particular, with respect to phase transitions 
and multifragmentation. We want discuss the question whether, starting from 
a transport description of a heavy ion collision, where in principle everything 
is known about the system, a thermodynamical picture of multifragmentation 
can be deduced. 

2 Determination of temperature 

In this work we make use of two different methods of determining temperature: 
In the first, we determine local temperatures by fitting the local momentum dis- 
tributions (obtained in our case from relativistic transport calculations) to co- 
variant Fermi-Dirac distributions at finite temperature in the local rest frame 
[0]. Non-equilibrium effects are taken into account by allowing a parametriza- 
tion of the momentum space distribution in terms of two thermalized Fermi 
spheres [|l| (or covariantly by ellipsoids). With this method we obtain a local 
microscopic temperature, Tioc- 

In the second method we follow the experimental method of fitting frag- 
ment energy spectra. These are generated in our calculations by applying 
a phase space coalescence algorithm |I[] to the final stages of the transport 
calculations. As in experimental analyses these spectra are interpreted in a 
Siemens-Rasmussen |^ or blast model 0, which assumes a thermalized freeze- 
out configuration of nucleons and fragments with a collective radial flow profile 
and a unique temperature |^, ^]. In this model the fragment spectra are given 
by 
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where n{(3) is the flow profile, 7 = \/l — (3"^ and where T is the global tempera- 
ture {a = 7/3p/T). The flow profile is obtained from the simulation. Then the 
remaining parameter is the temperature T which is fitted to experimental, resp. 
generated fragment spectra. In experimental analyses a global temperature is 
assumed, which characterizes the shapes of all fragment spectra. This is not 
obvious and should be clarified in the analyses of our transport calculations. 



3 Analysis of spectator matter (semi— central collisions) 

We first discuss the thermo dynamical properties of spectator matter in semi- 
central Au on Au reactions at intermediate energies. Tfiis reactions fias been 
studied extensively by tlie ALADIN collaboration 0. We determine the spec- 
tator temperature from fits to local momentum distributions (for more details 
see Ref. [^). Fig. 1 shows the time evolution of the temperature in the spec- 
tator (left side) for different beam energies. When the spectators are clearly 
developed in the transport calculations after about 40 fm/c, their tempera- 
tures approach a rather constant value of about T ^ 5 MeV which remains 
fairly stable up to about 80 fm/c, and furthermore is rather independent on 
the incident energy considered. 
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Fig. 1. Left: Temperature evolution in the spectator in semi-central Au 
on Au reactions at different beam energies indicated in the figure. Right: 
Density-pressure trajectories for the spectator matter for the same reaction at 
600 AMeV. The solid and dashed curves represent londitudinal and transverse 
pressure, respectively. The squares and circles are the values at different times 
starting from t = 35 fm/c in steps of 5 fm/c. The dotted curves are the nuclear 
matter isothermal equation of state for T = 5 and 7 (lower and upper curve, 
respectively). 

These results are in good agreement with experiments of the ALADIN collab- 
oration, which from measurements with different "thermometers" determine 
the same value of T Ri 5 MeV, depending only moderately on the beam energy 
of the reaction. We also generate pressure-density trajectories for the specta- 
tor matter as a function of time (right side of Fig. 1). Dynamical instabilities 
should arise when the pressure increases with decreasing density indicating 
a negative effective compressibility, which occurs here at t > 50 fm/c. The 
system at this stage therefore enters an instability region and should break up 
into fragments. Comparing to the nuclear matter isothermal equation of state 
for temperatures of T = 5 and 7 MeV corresponding to the range of spec- 



tator temperatures in Fig. 1 one sees that the thermodynamical conditions, 
as determined here, are close to but not identical to those of equilibrated nu- 
clear matter. Only at the final stages the spectator closely follows the nuclear 
matter behavior at temperatures of about T ^ 5 MeV. The densities at the 
instability condition are about 1/3 — 1/2 of saturation density. It thus appears 
that the spectator closely approaches a freeze-out configuration in thermal 
and chemical equilibrium. 

4 Analysis of the fireball (central collisions) 
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Fig. 2. Temperatures (left) and radial flow (right) obtained from blast model 
fits to fragment energy spectra as function of the beam energy. The theoretical 
results are shown for two mean field models (see text). The data are taken 
from I, ||. 

In central collisions the situation is rather different. If very central events are 
selected experimentally using charged particle multiplicities |^ or theoretically 
at polar angles near mid rapidity 0, there is no spectator matter. Rather 
one observes a hot dense fireball which expands isotropically as found by our 
calculations and also by other groups 0. Thus, assuming thermalization one 
can use Eq. (|I]) to extract the mean coUective radial fiow /?/=</?> and a 
slope temperature Tgiope from fits to the fragment energy spectra. Fig. 2 shows 
the energy dependence of these quantities as determined from our calculations 
and from experiments for central Au on Au collisions. Two parametrizations of 
the mean field (non-linear Walecka model and configuration dependent Dirac- 
Brueckner mean fields |jl[]) were used in the calculations to demonstrate the 
moderate dependence of of f]f and T on the mean field. As seen in Fig. 2 the 
experimental data for the radial fiow are reproduced very well. The comparison 
of the extracted slope temperatures Tgiope (left side of Fig. 2) is, however, only 
qualitative. 



It is of interest, to discuss the relation of these slope temperatures to the 
local temperatures Tioc determined from the momentum distribution of the 
calculation (we used a Maxwell-Boltzmann distribution here, in order to be 
consistent with eq. (|I]), but the difference is < 5 MeV in the final stages of the 
collision). It is also of interest to make the comparison separately for different 
fragment masses in order to determine whether a freeze-out szenario with a 
unique temperature is realistic. 
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Fig. 3. Temperatures (left) and radial flow (right) for fragments of different 
mass as determined from blast model fits. Also simultaneos fits to all frag- 
ments including nucleons (filled circle) and only to fragments with Aj > 2 
(filled triangle) are shown. The open diamond in the left figure is the lo- 
cal temperature at freeze-out obtained from a fit to local momentum space 
distributions (see text). 

The results of this analysis are shown in Fig. 3. Here we show the slope 
temperatures and radial flow velocities from blast model fits to the spectra of 
different fragment masses Af. The temperature increases and the radial flow 
decreases with increasing mass. In the coalescence picture such a behavior is 
reasonable, since a larger fragment has to be generated more inside the fireball, 
where the flow velocity is smaller and the temperature higher. In Fig. 3 we 
also show the result of a simultaneous blast model fit to all fragments and to 
fragments with mass Af > 2. Since the fragment multiplicities are roughly 
exponential and thus dominated by the nucleons the results for all fragments 
are close to those for Af = 1 alone. On the other hand the fit to the heavier 
fragments alone has lower radial flow and higher temperature and is the one 
compared in Fig. 2 to the corresponding experimental value. 

Also shown in Fig. 3 is the local temperature from the momentum distri- 
butions determined at about 35 fm/c. At this time the fireball in the calcula- 
tions approaches a freeze-out configuration (nucleon-nucleon collisions cease) 



in equilibrium (pressure isotropic). It is then a consistent check that the local 
temperature for this situation agrees approximately with the one determined 
from a blast model fit to the Af = 1 energy spectra. 

5 Conclusions 

We have studied the thermodynamical state of nuclear matter in heavy ion 
collisions by analyzing local phase space configurations and by analyzing frag- 
ment energy spectra. For the spectator a consistency with temperatures and 
breakup conditions with results of the ALADIN collaboration was found. For 
the participant matter we have applied in addition a blast model analysis to 
fragment spectra generated in the coalescence model. We see that the slope 
temperatures in such a description do not yield a unique value for all fragments. 
This does not favour the picture of a freeze-out configuration in thermody- 
namical equilibrium. Rather it appears that fragment emission is a dynamical 
process which occurs during a longer stage of the heavy ion reaction. 
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